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Flow  induced  by  a body  moving  near  a plane  wall  is  analysed  on 
the  assumption  that  the  normal  distance  from  the  wall  of  every  point 
of  the  body  is  small  compared  to  the  body  length.  The  flow  is 
irrotational  except  for  the  vortex  sheet  representing  the  wake.  The 
gap-flow  problem  in  the  case  of  unsteady  motion  is  reduced  to  a non- 
linear first-order  ordinary  differential  equation  in  the  time  variable. 
Problems  solved  include  airfoil  starting  flows  and  their  transient 
wakes,  and  flat  plates  falling  toward  the  ground. 
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1 . Introduction 

"'S>* Fluid -dynamic  problems  involving  bodies  moving  close  to  walls 
are  of  interest  in  many  different  contexts,  and  there  is  a considerable 
literature  dealing  with  such  problems.  The  present  paper  concerns 
itself  with  phenomena  which  can  be  treated  on  an  inviscid-fluid  basis, 
and  hence  has  little  connection  with  the  important  branch  of  that 
literature  dealing  with  low- Reynold s-number  wall  effects. 

The  classical  inviscid  problem  is  of  course  the  aerodynamics  of 
wings  near  the  ground,  and  indeed  the  present  work  has  its  most, 
obvious  application  in  that  area.  The  regime  of  interest  in  the 
present  paper  is,  however,  that  of  very  small  clearances  to  the 
ground.  The  early  work  on  aerodynamic  ground  effect  (e.  g.  as 
surveyed  by  Pistolesi,  1937),  views  all  effects  of  the  ground  as  small 
perturbations  to  the  infinite -fluid  flow  about  the  wing.  The  formal 
requirement  for  this  to  be  valid  is  that  the  clearance  between  body 
and  ground  be  large  compared  to  all  length  scales  of  the  body. 

An  intermediate  regime  in  which  the  flow  is  qualitatively  similar 
to  that  in  an  infinite  fluid,  but  where  ground  effects  are  "0(1)" 
perturbations,  is  when  the  clearance  is  comparable  to  some  body 
dimension.  For  example,  (e.  g.  Bagley,  1961,  Tuck  and  Newman, 
1974)  the  problem  of  a thin  airfoil  in  steady  ground  effect  at 
clearances  comparable  to  the  chord  (thus  large  compared  to  the 
thickness  and/or  camber),  leads  to  a singular  integral  equation  for 
the  bound  vorticity,  whose  kernel  reduces  to  the  classical  lifting - 
surface  kernel  for  an  unbounded  fluid,  as  the  clearance/chord  ratio 


tends  to  infinity.  On  the  other  hand,  if  one  lets  the  clearance /chord 
ratio  tend  to  zero  in  such  a theory,  the  integral  equation  "collapses", 
and  one  obtains  an  almost-trivial  explicit  result  for  the  loading, in 
terms  of  the  local  distance  from  the  wall. 

The  small-gap  regime  is  defined,  formally,  as  that  in  which 

^-7'  6i  /, 

the  clearance  is  small  compared  to  the  horizontal  length  scale/  For 

t 

non-thin  bodies  (e.  g.  automobiles), a general  approach  to  this  class 
of  problem  is  outlined  by  Tack  (1975);  however,  as  with  any  bluff- 
body  flow,  little  progress  can  be  made  without  introducing  empirical 
assumptions  regarding  the  wake.  We  assume  here  that,  in  addition, 
the  body  is  thin,  i.&  that  not  only  its  lower  surface,  but  also  its  upper 
surface,  is  close  to  the  wall. 

Some  literature  does  exist  on  the  small -gap  problem.  For 
example.  Strand,  Royce  and  Fujita  (1962)  noted  the  "hydraulic"  or 
channel-flow  character  of  the  tightly-constrained  flow  between  the 
body  and  the  wall,  and  Widnall  and  Barrows  (1970)  provided  a 
complete  asymptotic  solution  for  the  steady -flow  case,  assuming  in 
addition  that  the  thickness  and  camber  (and  "angle  of  attack  times 
chord")  are  small  compared  to  the  clearance.  The  present  paper  can 
be  considered  as  an  extension  of  the  work  of  Widnall  and  Barrows  to 
include  "non-linear"  (thickness  etc.  comparable  to  clearance)  and 


unsteady  effects.  Unsteadiness  due  to  ground  irregularity  was  included 
in  the  linear  theorv  bv  Barrows  and  Widnall  <1970i. 
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Although  the  potential  applications  are  to  very  practical  problems, 
such  as  tracked  ground  vehicles  (Barrows,  1971).and  ship  manoeuvring 
in  shallow  water  near  banks  (Norrbin,  1974;  Tuck,  1978),  in  the  present 
paper  we  give  only  simplified  illustrations,  involving  two-dimensional 
flow,  and  specialized  geometry  and  motions.  The  small-clearance 
assumption  implies  that  a two-dimensional  flow  becomes  one-dimen- 
sional in  the  gap,  and  is  thus  describable  by  the  one-dimensional 
continuity  equation.  For  a given  body  under-surface,  this  equation  is 
a second-order  ordinary  differential  equation  for  the  velocity  potential, 
as  a function  of  the  horizontal  coordinate  x along  the  wall,  which  can 
be  solved  explicitly.  Two  boundary  conditions  are  needed,  and  these 
must  come  from  matching  with  the  outer  flow  passing  over  the  top  of 
the  body.  In  the  Appendix,  we  show  that  the  proper  conditions  are 
continuity  of  velocity  potential  at  the  leading  edge,  and  (for  bodies 
without  stern  appendages)  of  pressure  at  the  trailing  edge. 

When  these  conditions  are  applied  in  unsteady  flow,  the  result  is 
a non-linear  ordinary  differential  equation  of  the  first  order  with  re- 
spect to  time,  to  determine  the  net  flux  exiting  from  the  gap  at  the 
stern.  In  the  special  case  of  steady  flow,  the  appropriate  solution  of 
this  equation  is  vanishing  flux,  in  a frame  of  reference  fixed  in  the 
fluid  at  infinity.  That  is,  in  a frame  of  reference  fixed  in  the  body, 
the  trailing-edge  velocity  is  equal  to  the  free-stream  velocity.  The 
present  theory  then  reduces  to  a non-linear  equivalent  of  that  of 
Widnall  and  Barrows  (1970).  Some  non-linear  consequences  for  the 
steady  lift  force  are  discussed  in  $4. 
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In  §5  we  analyse  the  transient  establishment  of  the  steady  flow 
of  §4,  starting  from  rest.  The  differential  equation  for  the  gap  flux 
at  the  6tern  can  be  solved  explicitly,  for  the  case  of  an  impulsive 
start  at  constant  speed.  In  this  simple  case,  we  can  also  solve 
for  the  velocity  and  thickness  of  the  vortex  wake  left  behind  the  body, 
and,  for  example,  predict  a shock-like  rolling  up  of  the  wake  behind 
an  airfoil  at  a negative  angle  of  attack. 

The  problem  solved  in  §6  is  free  fall  under  its  own  weight  of  a 
flat  plate  toward  a plane  wall,  combined  with  a horizontal  motion  at 
constant  speed  U.  For  U = 0 , a symmetric  solution  was  presented 
by  Yih  (1974),  in  which  the  plate  never  actually  hits  the  wall,  but 
ultimately  approaches  it  with  an  exponential  decay  of  height  with  time. 
We  find  a similar  result,  but  only  for  U above  a certain  critical 
value.  For  lower  values  of  U,  the  approach  to  the  wall  is  still 
gentle,  but  there  is  impact  at  a finite  time,  with  a cubic  decay  of  height 
with  time  near  impact.  It  should  be  noted,  however,  that  we  retain 
the  asymmetry  between  leading  and  trailing  edges  even  when  U is 
small,  and  Yih's  symmetric  solution  is  likely  to  be  more  realistic  for 
small  U.  We  also  assume  that  the  angle  of  attack  remains  zero  for 
all  time,  and  (at  least  when  U ± 0)  this  makes  the  application  to 
problems  such  as  stacking  of  glass  plates  (Yih,  1974)  or  the  "sliding 
of  sheets  of  paper"  illustration  in  G.  I.  Taylor's  movie  (1967)  not 
yet  complete.  It  would  appear,  however,  that  there  is  no  need  to 
invoke  viscous  effects  to  explain  a number  of  air-lubrication 
phenomena. 


s' 
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Indeed,  in  many  respects  the  present  paper  merely  opens  the 
door  to  further  investigation  of  interesting  small-gap  phenomena. 

The  potential  applications  to  ship  manoeuvring  are  particularly 
important.  Although  we  have  phrased  the  theory  here  in  terms  of 
wall  effects,  there  is  also  no  difficulty  in  extending  it  to  the  more 
general  problem  of  two-body  interactions,  e.  g.  ships  in  passing 
manoeuvres,  or  biplane  wings. 

Extensions  to  three-dimensional  flow  can  take  a number  of  forms. 
Perhaps  of  greatest  practical  interest  is  the  case  where  a ship  is 
moving  in  shallow  water  with  a very  small  (compared  to  its  beam' 
clearance  to  the  bottom,  and  a somewhat  larger  clearance  (e.  g. 
comparable  to  its  beam)  to  a side  wall  or  another  ship.  In  that  case, 
the  one -dimensional  assumption  remains  valid  between  ship  and  wall, 
but  there  is  possible  leakage  of  fluid  under  the  ship,  i.  e.  out  of  the 
"channel".  This  geometry  is  similar  to  that  considered  by  Barrows 
(1971)  for  a ram  wing  in  a "trough". 
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2.  Problem  Formulation 

We  assume  two-dimensional  irrotational  flow  of  an  incompressible 
fluid,  generated  by  movement  of  a thin  airfoil -like  body  as  sketched 
in  Figure  1,  with  upper  surface  y = f , and  lower  surface  y = f , 
i.  e.  occupying  the  region 

fL  (x,  t)  < y < fu(x,  t)  , xR(t)  < x < Xg(t)  . (2.1) 

The  leading  edge  or  bow  is  at  x = xR  , and  the  trailing  edge  or  stern 
at  x = Xg  , The  body  has  length 

~ ~ Xg  , (2.  2) 


which  can  in  the  most  general  case  depend  on  time,  but  is  normally 
considered  to  be  constant. 

The  body  is  everywhere  close  to  a plane  boundary  at  y = 0 , i.  e. 


fL’  fU  * °«'> 


(2.  3) 


where  « is  a small  parameter.  This  also  requires  the  body  to  be 

thin  , as  in  the  thin-airfoil  theory.  However,  we  do  not  assume  that 

its  thickness  fy  - f ^ is  small  compared  to  the  wall  separation,  e.  g. 

to  the  minimum  value  of  f . 

L 

Our  task  is  to  solve  Laplace's  equation 


<f>  4 0 

**  ™yy 


0 


(2.  4) 
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for  the  velocity  potential  ♦ (x.y.t)  . subject  to  suitable  boundary 
conditions.  At  infinity  we  have  a state  of  rest.  i.  e. 

0 . V p — 0 as  x,  y «o . (2.  «») 

The  wall  y 0 is  inipe rmeable, i.  e. 

* lx.  0.  t)  0 . (2.  M 

V 


The  boundary  conditions  on  the  moving  body  surface  are 


0 f + 
y 


(2.  7) 


where  we  use  a dot  for  an(^  A prime  for  • Equation  i2.  7) 

applies  both  with  f f^  and  f = f^  . 

Finally  we  must  use  an  appropriate  wake  and  Kutta  condition.  The 
trailing  edge  \ is  assumed  to  shed  vortices  which  remain  behind 

the  body  in  a vortex  sheet,  with  equation 


y - fyyiX.t)  . 


(2.  81 


the  function  f^(x.t)  being  unknown.  The  kinematic  boundary  condition 


across  this  surface  is  given  again  by  (2.7),  with  f = f^  . 


The 


dynamic  condition  is  continuity  of  pressure  across  the  sheet.  Thus  if 
p(x.  y,  t)  is  the  excess  of  pressure  over  the  value  at  infinity,  then  from 


f 

i 
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Bernoulli's  equation 


P = -P(0  + £ ^ + I <t>  ) 


(2.  9) 


everywhere  in  the  fluid.  The  dynamic  boundary  condition  on  the  wake 


p(x,  + 0,  t)  = p(x,  - 0,  t)  . 


(2.  10) 


Equation  (2.  10)  holds  in  the  wake  for  x > xg(t)  and  also  holds  at  the 
trailing  edge  x = *g(t)  » where  it  defines  the  Kutta  condition,  and 
ultimately  determines  the  circulation  about  the  body. 
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One -Dimensional  Theory 

We  now  assume  that  in  the  limit  as  < - 0 . $ 0(«  ) except 

in  the  near-wall  gap  region  G defined  by 

G:  0 < y < f(x,  t)  , x > Xg(t)  , (1.1) 


whe  re 


f *L(x.  t)  , 


f(x.t) 


xB(t)  < X < xs(t) 
X > x (t)  . 


( 1.  £) 


In  the  region  G between  the  body  (or  its  wake!  and  the  wall,  there  is 
a flow  of  a magnitude  which  does  not  tend  to  a>ro  as  * — 0 . if  We 
are  interested  only  in  leading-order  estimates  of  forces,  etc., we  can 
therefore  concentrate  attention  on  the  flow  in  region  G , assuming 
m effect  that  $ = 0 when  we  are  outside  G . A justification  for  this 
assumption  is  provided  in  the  Appendix  by  matching  techniques  (see  also 
Widnall  and  Barrows,  1970). 

In  tact,  the  flow  in  G is  classical  one -dimensional  or  channel 
flow,  in  which  the  fluid  moves  predominantly  in  the  x direction, 
parallel  to  the  wall.  For  example  it  is  consistent  with  this  approximation 
to  use  the  Taylor  expansion 

= tf>(x,0,t)  - jy2  (x.  0,t)  + (3.1) 


auaiMHIfil 


1 1 


which  guarantees  satisfaction  of  (2.  4)  and  (2.  6),  and  satisfies  (2.  7) 
if 


f + (f  / )'  = 0 . 


(3.  4) 


Alternatively,  one  may  recognize  (3.  4)  as  the  ordinary  one -dimensional 
continuity  equation,  expressing  conservation  of  mass  in  a flow  with 
dominant  velocity  <b'  (x,  0 , t ) , in  a channel  of  width  f which  varies 
both  in  space  and  time.  From  now  on,  we  write  </>(x,  t)  for  <t>(x,  0,  t)  . 

One  integration  of  (3.  4)  gives 


f(x,  t)4>'  (x,  t)  = q(t ) + 


1 


xs(t) 


f(4,t)d£  , 


(3.  5) 


where 


q(t)  = f(xs(t),  t)0,(Xg(t),t) 


(3.  6) 


is  the  (so-far  unknown)  net  flux  through  the  gap  at  the  stern.  Our 
primary  task  is  to  determine  this  quantity  q(t)  . A further  integration 
gives  the  potential  <p  itself.  In  performing  this  integration  step,  we 
match  with  the  exterior  flow  0 = 0 outside  G , by  applying  the 
boundary  condition 


[1 


0(xB(t),t)  = 0 , 


(3.  7) 
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at  the  leading  edge  x --  Xg(t)  . This  i8  justified  formally 
Appendix.  Thus 


Mx,t)  = q(t)A(x,  t)  + B(x,  t)  , 


where 


A 


/ 


1 

7 


i. 


e. 


A(x, t) 


■c 


(t) 


di 

Ct) 


and 


B 


fdl  , 


e. 


- r -ai_  c ,s(t) ... 

J f(4.t)  J . 


B(x,  t) 


t)d| 


in  the 


(3.  8) 


(3.  9) 


(3.  10) 
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For  a given  body  lower  surface  f = f^  , (3.  8)  determines  the 
flow  in  that  part  of  G beneath  the  body,  providing  we  can  find  q(t)  . 

In  order  to  find  q(t)  , we  must  make  use  of  the  Kutta  condition  (2.  10) 

at  the  trailing  edge  x = x (t)  . Now  in  G , we  can  simplify  the 
Bernoulli  equation  (2.9)  since  <t> y = 0(c  ) , and  hence  the  pressure 
in  G is  given  by 

p(x.t)  = -p[*(x,t>  + |V(x.t»2]  + 0(*2).  (3.11) 

Since  4>  = 0 outside  G , and  hence  p s.  0 (actually  p = 0(« ))  , 
continuity  of  p across  the  wake  simply  requires  p to  vanish  there, 
i.  e. 

4>(x,t)  + ](4>'(x,t))2  - 0 , x > xg(t)  . (3.  12) 

Equation  (1.  12)  is  a first-order  non-linear  partial  differential  equation 
to  determine  the  flow  in  the  wake  region,  equivalent  to  the  homogeneous 
Euler  equation 

4 u ^ = 0 (3.  13) 

9t  9x 

for  the  velocity  u Once  (3.  12)  or  (3.  13)  is  solved,  we  have 

another  first-order  equation  (3.  *1)  to  solve  for  the  unknown  wake 

boundary  y f (x.t)  . Such  a computation  is  presented  in 

w 

However,  it  in  not  necessary  to  go  this  far,  if  we  are  primarily 
interested  only  in  the  flow  beneath  the  body. 
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So  long  as  the  gap  and  wake  regions  arc  essentially  continuous 
with  each  other  (i.  e.  not  interrupted  by  an  edge  region  where  ^ is 
significant),  equation  (3.  12)  can  also  be  applied  "at"  x = Xg(t)  , 
and  is  the  required  Kutta  condition  to  determine  the  trailing -edge  flux 
q(t)  . In  the  Appendix  we  show  that,  if  the  under-surface  at  the 
trailing -edge  is  "shaped",  or  if  there  is  an  appendage  such  as  a rudder 
at  a finite  angle  of  attack,  (3.  12)  is  replaced  by 


* + yz 


(3.  14) 


where  o is  the  contraction  coefficient  of  the  "nozzle"  formed  between 
the  trailing  edge  and  the  ground  plane.  If  o = 1 , as  we  assume 

henceforth,  (3.14)  reduces  to  (3.  12)  at  x = Xg- 

On  substitution  of  the  representation  (3.  8)  for  the  velocity  potential 
<p  into  the  Kutta  condition  (3.  12)  at  x = Xg(t)  , we  obtain  the  non- 
linear first-order  ordinary  differential  equation 

2 

A(xc.t)q  4 A (xc,  t )q  + B(x<;,t)  + j- ^ - 0.  ( 3.13) 

S S b 2f  (Xg.t) 

• • 

Note  that  the  coefficients  A,  B are  obtained  by  differentiating  the 
corresponding  expressions  (3.  9),  (3.  10)  with  respect  to  time,  and 
then  setting  x = xg(t)  ; in  particular,  A(xg,t)  i A(Xg,t). 

The  differential  equation  (3.  15)  must  be  solved  in  any  special  case, 
subject  to  a suitable  initial  condition  on  the  flux  q(t)  at  t = 0 . 


m 


The  above  applies  for  an  arbitrarily -deforming  body  surface. 
We  are  most  interested  in  the  case  of  rigid  bodies  with  3 degress 
of  freedom,  measured  by  the  co-ordinates  (xQ(t),  y Q(t ) ) of  the  mid- 
point of  the  body,  and  the  (bow-up)  pitching  angle  0(t)  about  that 
point.  If  (X,  Y)  are  co-ordinates  fixed  in  the  body.aligned  with  the 
(x,  y)  axes  at  6 = 0,  and  if 


Y = fQ(X)  , | X | < f . 


(3.  16) 


is  the  equation  of  the  lower  surface  of  the  body  in  these  co-ordinates, 
then  when  motion  occurs,  the  lower  surface  is  given  by  y = f(x,t)  . 
whe  re 


f(x,  t)  - yQ(t)  - 0(t)(x  - xQ(t))  + fQ(x  - xQ(t))  . 


[3.  17) 


On  substitution  into  (3.  9),  (3.  10),  we  have 


A(x,  t) 


(3.  18) 


B(x,t) 


Vnri  fl?<. + eiow  -6) 


-\6iZ  - x f (f) 
o o 


+ + xQf  (g)l 


(3.  19) 
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4.  Steady  Flow 

An  important  special  case  is  when  xq  = -Ut  . yQ  = 0 = 0 . 
That  is,  the  rigid  body  with  lower  surface  y = f^(x  + is  m°ving 
with  constant  velocity  U in  the  negative  x-direction.  Although  it  is 
somewhat  easier  to  solve  this  problem  directly  in  the  moving  (X,  Y) 
co-ordinate  system,  in  which  the  body  appears  fixed  in  a uniform 
stream  of  magnitude  II  in  the  +x  direction,  we  provide  the  solution 
here  by  specialization  of  the  general  unsteady  results. 

We  observe  that 


A 


I 


x + Ut 


dX 

foon 


(4.  I) 


and 


B Ufo(/  )A  - (x  + Ut  + f )U  , 


(-1.  2) 


and  find  that  B - 0 at  the  stern.  Hence  the  differential  equation  (3.  15) 
possesses  the  solution  q =•  0 , which  simply  means  that  free-stream 
conditions  prevail  at  the  stern.  Equation  (3.  8)  then  shows  that  (4.  2) 
is  in  fact  an  expression  for  the  velocity  potential  0(x, t)  everywhere 
in  the  gap.  from  which  we  can  compute  (using  (3.  11))  the  pressure  in 
the  gap, 


fo('> 

fJoT+ut) 


l 


(4.  3) 
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Since  we  have  assumed  that  p - 0(«  ) on  the  upper  surface  of  the 
body,  the  net  upward  force  per  unit  span  is 


Fl  = pU Z£(l  - \ ) + 0(c  ) , 


(4.  4) 


whe  re 


X(x)  = fQ(x)/fo(f) 


(4.  5) 


is  the  gap  thickness,  scaled  with  respect  to  r.hat  at  the  trailing  edge, 
and  a bar  denotes  a mean  value  over  the  length  of  the  body,  i.  e. 


1 

2 i 


dx 

(Mx))2 


(4.  6) 


The  result  (4.  4)  is  equivalent  to  one  obtained  by  Strand,  Royce  and 
F ujita  (1962). 

It  is  instructive  to  compare  the  above  lift  force  with  that  which 
would  be  obtained  by  use  of  the  Kutta-Joukowski  theorem,  which 
applies  only  in  unbounded  media.  Thus  (since  <f>  is  continuous 
elsewhere),  the  net  circulation  T around  the  body  is  generated 
entirely  by  the  jump  in  the  velocity  potential  <t> , from  its  value 
^(Xg,t)  immediately  below  the  trailing  edge,  to  its  zero  value  immediately 


I 


above  it.  Hence 
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r $ixs,t) 

= -2U'  + ufo«'>  r rn 

-I  ° 

= -2U/(1  - X71"  ) . 7) 

Thus  the  Kutta  - Joukowski  lift  is 


- (>  U r 

2PU2|(1  - X*1"  ) . 


(•».  8) 


In  cases  where  the  body's  thickness  or  camber  is  small  compared 
to  the  wall  separation,  the  Kutta -Joukowski  theorem  becomes 
asymptotically  valid,  and  F'  •*  F.  . since  in  that  limit  (f  (x)  = f (/)) 

K Lj  O O 


2(1 


- 2(Y  - l ) . 


(-».  11 ) 


That  is. 

- 2pU2<  (5T  - 1)  , (.1.  io) 

which  depends  only  on  the  net  area  2 1 (1C  - 1)  lying  between  the 
lower  surface  of  the  body  and  a horizontal  line  through  its  trailing 
edge.  For  example,  if  the  body  is  a flat  plate  at  an  angle  of  attack  a . 


then 
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F L = |pU2.  Zi.a.  ( j^jj)  , 

and  is  inversely  proportional  to  the  clearance /chord  ratio. 
The  result  (4.  11)  was  found  by  Widnail  and  Barrows  (1970). 


(4. 
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5.  Impulsively -Started  Body  at  Zero  Angle  of  Attack 

We  assume  again  that  6 = 0,  yQ  = 0 , and  hence  the  gap 
thickness  is 

f(x,t)  = fQ(x  - xQ(t))  . (5.  1) 

The  differential  equation  (3.  15)  becomes 

Aoq  - - W +-4=0.  (5.2) 

o 2f 

o 

where 

Ao  = I = u X_1  /£0(i)  * <5- 3) 

and 


f = 
o 


(5.  4) 


For  a general  pattern  x (t)  of  horizontal  motion,  we  must  solve 

o 

(5.  2)  numerically.  However,  for  the  special  case 


• • 
x 


o 


- U6(t)  , 


(5.  5) 


in  which  there  is  a sudden  commencement  of  steady  motion  toward  the 
left  at  speed  U , we  have  to  solve 
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V - f o >r?r 

o 


= 0 , 


(5.  6) 


subject  to  the  initial  condition 
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q(o)  = U(-~  - f ) , 

o ° 


(5.  7) 


which  follows  by  integrating  (5.  2)  from  t = 0-  to  t = 0+  . The 
solution  is  given  by 


2Uf 


o _ t/t  . 

- - fie  o - 1 , 


(5.  8) 


where 


H-  = (1  + X-1  ) / (1  -X.-1  ) # 


(5.  9) 


and 


f A /U 
o o 


= 2i  X." 1 /U 


fo(i) 


u 


I i 


dl 


.(&) 


(5.  10) 


As  t oe  , q — 0 , and  the  steady  flow  of  §4  is  recovered.  The 


! 


, 


time  scale  for  development  of  the  steady  flow  is  measured  by  t 
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e.  g.  if  X.  1 = i ( M-  = 3)  , then  the  flux  has  fallen  to  28%  of  its 
starting  value  when  t = t . It  should  be  noted  that  during  any  time 
interval  t, the  body  moves  a fraction  X. ’ * of  its  length.  However, 
if  it  happens  that  A.  * - 1 (p  - oc)  , the  steady  flow  is  established 
immediately,  since  q(o)  = o implies  q = o for  all  t > o . This  applies 
to  some  non-trivial  flows,  as  well  as  to  the  trivial  zero-disturbance 
case  X.(x)  = 1 , of  a flat  bottom  without  angle  of  attack. 

In  the  present  case  it  is  straightforward  to  solve  (3.  13),  (3.  4) 
for  the  wake  velocity  $' = u(x,t)  and  wake  width  f = f^(x,t)  , the 
problem  being  essentially  identical  to  that  for  a compressible  gas  at 
constant  pressure,  as  described  for  example  by  Ames  (1965,  p.  51), 
with  fyj  playing  a role  analogous  to  the  gas  density  p . Thus  the 
general  solution  of  (3.  13)  is  given  by  the  transcendental  equation 

x = ut  + G(u)  , (5.  1 1 ) 

for  some  function  G(u)  . Applying  the  boundary  condition 

u = 3151  at  x = - Ut  , (5.  12} 

o 

where  q(t)  is  given  by  (5.  8),  we  find 

G(u)  = - tp  (U  + u)  log  |2--^  | . (5.13) 

The  corresponding  general  solution  of  the  linear  equation  (3.  4) 


for  f is 


for  some  function  H(u),  determined  from  the  boundary  condition 


f 


W 


f at  x 
o 


- Ut 


(5.  15) 


to  be 


H(u) 


fQ  [G'(u) 


G(u)  i 
U + u J 


= 2f  Ut 
o o 


U + u 
u(2U  + u) 


(5.  16) 


In  Figure  2,  we  plot  u and  f^.  against  x,  for  various  values  of 
t,  with  p = 3 , i.  e.  X~  ^ \ . This  is  a case  with  effectively- 

positive  angle  of  attack,  and  hence  acceleration  of  the  fluid  in  the  gap. 
There  is  positive  stern  flux  for  all  time,  and  hence  the  vortices  shed 
at  the  trailing  edge  move  to  the  right,  and  the  wake  continually 
lengthens  and  thins. 

A very  different  picture  emerges  when  p is  negative,  or  X * > 1, 
i.  e.  when  there  is  an  effective  negative  angle  of  attack.  Now  there  is 
deceleration  in  the  gap,  and  negative  flux  at  the  stern.  Hence  the 
shed  vortices  move  to  the  left,  and  the  wake  rolls  up  on  itself.  This 
is  displayed  in  the  mathematic sfby  the  possibility  of  a singularity  in 
(5.  14)  if  t = - G'  , which  can  happen  only  if  p < 0 . The 
phenomenon  is  similiar  to  shock  waves  in  gas  dynamics. 
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6.  Falling  Flat  Plate  in  Steady  Horizontal  Motion 


We  now  assume  f (x)  = 0,6  = 0 , and  x 
o o 


-Ut  , but  allow 


an  arbitrary  vertical  movement  y = yQ(t)  • The  velocity  potential 
(3.  8)  becomes 


u(x  + Ut  + I ) + j P(x  + Ut  + 1 )(x  + Ut  - 31  ) 


(6.  1) 


where 


u = q/yQ  = /(-Ut  + i , t) 


(6.  2) 


is  the  gap  velocity  at  the  stern,  and 


p = - y /y 

^ 7o  7o 


(6.  3) 


The  differential  equation  (3.  15)  becomes 


• 2 • 1 2 
21  u - 21  cp  + Uu  + j u = 0 . 


(6.  4) 


For  any  given  time  history  of  vertical  movement  y (t)  , we  must 
evaluate  the  quantity  p(t),  and  then  solve  the  differential  equation  (6.  4) 
for  u(t)  • 

The  pressure  on  the  plate's  under  surface  is  given  by 


- £ = u (x  + Ut  + 1 ) + j p (x  + Ut  + 1 )(x  + Ut  - 31 ) + U[  u + (3(x  + Ut  - 1 )] 
P ^ 


+ 2 [ u + (3(X  + Ut  - 1 )J  2 , 


(6.  5) 
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and  vanishes  as  required  at  x + Ut  = i , if  (6.  4)  is  satisfied.  The 
net  upward  force  is 


pirziu  - 3i2p  + 2U(u  -pf  ) + (u  -pi  )£  + ^(P+P^)]  (6.  6) 


2 1.2, 

3 J 


and  the  (bow -down)  moment  is 


M - -f  pi  Vu  - ip  + (U  + u)p  - ip?1 


(6.  7) 


For  example,  consider  a freely -falling  plate  of  mass  m per 
unit  span,  constrained  against  rotation.  Then  y Q(t ) itself  is 
determined  by  solving  the  equation  of  motion 


Fl  = m(yQ  + g)  , 


(6.  8) 


i.  e. 


21  u 


- | i 2p  + 2 Uu  + u2  -2(U  + u)pi  + | i 2p2 


■T-  (?o  + *>  • 


(6.  9) 


Equations  (6.  3),  (6.  4)  and  (6.  9)  are  a set  of  3 coupled  non-linear  first- 
order  differential  equations  to  solve  for  u,  p and  yQ  . 


If  we  define  t = t as  the  instant  at  which  the  plate  hits  the  wall, 
o 


so  that  motion  takes  plate  for  t < t , then  so  long  as  t remains 


finite,  the  final  stage  of  the  fall  appears  to  be  described  by  the  estimates 


' II 


n 


u 


li 

t -t 


U - 


+ 11 


l 

~U 


(t-»o)  4 0(t-t  ) . 


(6.  ! 0) 


- t 


t -t 


Ufc(t-t  ) 
o 

■U 


+ °‘t"*or  * 


(6.  I 1 ) 


a nd 


* u2(t-t  r 

(t-t  ) [l  -t  y?—  1 0(t-t  ) | , 

8/  ° 


(6.  k>) 


l°r  some  constant  c . this  constrasts  with  the  exponential  approach 
to  the  wall  found  by  Yih  ( 197-1)  at  U 0 , by  making  the  empirical 
assumption  that  p 0 at  hot h edges  of  the  plate.  This  means  that 
vortex  sheets  must  spring  from  both  edges,  there  being  of  course  no 
distinction  between  "leading  " and  "trailing"  edges  if  U 0.  Although 
the  result  is  different,  the  conclusion  that  the  approach  to  the  wall 
is  relatively  'gentle"  still  holds.'  However,  the  requirement  of  no 
rotation  is  rather  restrictive,  since  the  limiting  bow  down  moment  is 

s r‘ 

M -T-  — t’ -t  Oil).  I (,.  l t) 

(t-tor 


In  tact  it  is  possible  to  have  an  exponential  approach  to  the  wall, 

for  sufficiently  - large  horizontal  speed  U . Thus,  the  system  (fe.  <). 

(6.  1),  (6.  ")  allows  a solution  with  u •*  0 , fl  — ft  , and 

’ ’ o 


-P  t 

r o 


y - 


i*  l% 


“*  0 , t ■*  or  , 


(b.  1») 
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whi're  is  .1  positive  constant  determined  from  tin*  quadratic 

equation 


*12  2 

-2  Ui  p t - t p 
o \ o 


i . e . p 


<6.  l^) 


(6.  U>) 


This  solution  possesses  a finite  bow-up  limiting  moment  M. 

The  exponentially-decaying  solution  can  exist  only  if  the  square 
root  in  (6.  16)  is  real.  i.  e.  if 


U > - 

V 3 of 


(6.  17) 


The  value  of  p^  at  this  limiting  speed  U is  one  half  of  the  value 
found  by  Yih  at  U 0 . and  as  U *> , the  dominant  (smaller  p ) 
decaying  term  has  p **  0 . The  present  theory  must  break  down 
for  sufficiently  small  U , and  Yih's  theory  appears  a suitable 
extrapolation  to  II  0 , with  the  advantage  of  flow  symmetry  about 
the  centre  of  the  plate. 

The  above  system  was  solved  numerically  for  motion  starting 

from  rest  at  y h . We  use  a non-dimensional  version  of  t.he 
7o 


equations,  namely 


29 


- Pv, 


, y(o)  = l . 


* • 1 rT  12 
u = p - T Uu  - -r  u 


, u(o ) = 0 , 


P = 


1 1 2 2 2 1 2 
- 4-  Uu  + UP  - -u  + u|3  - - p - y my  p - m 
£ ‘i  .5  £ O 


] 


l 


, p(o)  = 0,  (6.  18) 


-3-2  myo 


where  the  scale  for  time  is  the  time  taken  to  fall  a distance  h 
vacuum,  i.  e.  T 


/ 2h 
v K 


in  a 
(6.  19) 


and  y is  scaled  with  h,  u and  U with  f / T , p with  1 /T  , and 

3 

m with  pi  /h  . 

It  appears  from  these  computations  that  whenever  the  condition 

(6.  17)  is  satisfied  (i.  e.  high  speed  U) , the  exponential -decay  solution 

holds,  with  the  smaller  of  the  two  values  of  po  in  (6.  16).  On  the  other 

hand,  if  (6.  17)  is  not  satisfied  (low  speed  U),  the  body  does  indeed 

hit  the  ground  at  a finite  time  t , and  t -*■  oo  as  U “*•  U . 

" o o o 

Figure  4 shows  a typical  set  of  numerical  results  at  (scaled)  m - 1,  for 
various  values  of  (scaled)  U . The  critical  scaled  speed  is  Uf  ( = n/b/  3 
1.  64  in  this  case.  For  example,  we  find  at  U = 1 that  the  body 
hits  the  ground  (with  zero  velocity  and  acceleration  as  predicted  by 
(6.  12))  at  t - 3.  3 , whereas  at  U = 2 there  is  an  exponential  approach 
as  in  (6.  14),  with  p^  - 0.  63,  as  predicted  by  (6.  16).  The  curve  with 
U 0 has  no  physical  significance  but  is  included  as  a limiting  case 
for  small  U . Figure  4 also  shows  (dashed)  the  parabolic  free  fall 
in  a vacuum,  equivalent  to  (scaled)  m ■*  « for  all  finite  U.  Fluid- 
dynamic  effects  always  postpone  the  impact  (sometimes  forever!  ) , 
and  also  make  the  impact  gentle  if  it  does  occur. 
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Appendix 

Edge  -Region  Matching  Problems 

The  inner  one -dimensional  flow  that  is  the  main  topic  of  the 
present  paper  is,  in  an  asymptotic  analysis  of  the  complete  flow, 
matched  with  an  outer  flow  above  the  body  via  separate  edge-region 
flows  near  the  bow  and  stern.  For  linearized  steady  flows,  the  details 
have  been  given  by  Widnall  and  Barrows  (1970),  and  the  present 
extension  to  non-linear  unsteady  flows  does  not  alter  the  fundamental 
character  of  the  matching  procedure. 

In  fact,  it  is  clear  on  stretching  that  the  flow  in  the  edge  regions 
is  quasi-steady,  in  a frame  of  reference  moving  with  the  velocity  of 
the  end  point.  This  is  because  spatial  rates  of  change  dominate 
temporal  rates  of  change  by  a factor  of  0(«  * ) in  the  edge  regions. 

For  example,  the  flow  at  the  bow  can  be  written  as 

0(x,  y,  t ) - *B(t)  - h(t)  uQ(t)  «(  - — ) (Al) 

where  #(X,  Y)  is  a canonical  non-dimensional  steady -flow  potential, 

as  sketched  in  Figure  5,  and  , u^  are  the  apparent  limiting  values 

of  <p  , , as  seen  by  an  inner  observer  as  x *•  x.,  4 o.  The  actual 

ox  n 

gap  size  at  the  bow  is  written  as 

hit)  = fL(xB(t),  t)  . (A2) 

and  this  is  scaled  to  a unit  value  for  the  non-dimensional  potential 
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#(X,  Y).  Notice  that  fc^XasX-^-oo,  o < Y < l , so  that, 

on  matching  with  the  inner  flow  as  x •*  x + o , 

B 


<t>(x,t)  ~ <(>B  + uB(x  - xQ)  . 


/?  2 
X + Y "*  ao , 

X - - oo,  o < Y < 1 , 


(A3) 


everywhere  other  than 


$ 


n 


log(Tie 


+ Y2  ) 


A4) 


(the  constant  "we"  being  obtained  by  conformal  mapping  solution  as  in 
Widnall  and  Barrows,  1970).  Thus,  as  JxZ  + y2  - oo  , 

hu„  hun  / 

* ((pB  ' log  TT}  ‘ T log /(x-x b)  + y • (A5) 

which  represents  an  apparent  sink  at  (xg  , o)  , as  seen  by  an  outer 
observer,  with  the  correct  strength  to  generate  the  net  flux  hu 

B 

entering  at  the  bow.  Since  the  outer  flow  is  everywhere  0(*  ),  the 
"constant"  term  in  (A5)  can  only  be  of  order  t at  most.  Since  h = 0(*  ) 
already,  this  means  that  #B  = 0(e  ) , and  hence  we  can  assume  as  our 
leading. edge  condition  on  the  inner  potential  0(x,  t)  that  0(x  , t)  = 0 , 

as  in  the  text.  This  condition  was  also  used  by  Barrows  and  Widnall 
(1970). 

Figure  5 is  drawn  for  leading  edges  without  0(1)  curvature. 
However,  curved  leading  edges  do  not  alter  the  qualitative  features 
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of  the  above  argument,  only  the  value  of  the  constant  replacing  "no" 
in  (A4).  In  particular,  the  conclusion  ■-  ()(*  ) still  applies  for 
curved  leading  edges,  so  long  as  the  curvature  remains  0(1)  only 
within  0(t  ) of  the  edge. 

The  flow  near  the  stern  is  not  like  that  in  Figure  5,  since  we 
assume  that  the  stern  flow  separates  smoothly  into  a wake,  rather 
than  negotiating  the  sharp  edge  to  become  source -like  at  infinity. 

In  fact,  if  the  trailing  under -surface  does  not  possess  0(1)  edge 
curvature,  the  stern  flow  is  extremely  simple,  consisting  of  nothing 
more  than  a uniform  stream  emerging  from  beneath  a parallel  plane, 
with  a constant  wake  height  equal  to  the  trailing -edge  height  f (xc.  t)  . 
I hus,  if  x - Xg  0(«  ) , the  stern  edge  flow  is  given  by 

0(x,y.t)  *s(t)  + us(t)(x  - xg(t))  + 0(«2)  ( A6) 


which  matches  the  inner  flow,  with  , ug  as  the  apparent  trailing- 
edge  values  of  respectively.  The  wake  condition  (2.  10)  for 

x > xg  reduces  to  (3.  12)  in  this  case,  since  p(x,  fw  + o.t)  and  <t>^  o , 
and  is  satisfied  with  0(c  ) error,  providing 


' USXS 


, 1 2 n 
+ 2 US  ‘ 0 * 


(A7) 


This  is  equivalent  to  the  Kutta  condition  for  the  inner  flow,  i.  e.  to 
(3.  12)  applied  at  x = Xg(t)  - 0 , since 


b 
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= ^-  0(Xg(t),t)  = ^(Xg(t),t)  + UgXg  . (A8) 

Somewhat  more  interesting  is  the  case  where  the  lowe resurface 
geometry  near  the  stern  is  curved  to  0(1)  slopes,  as  sketched  in 
Figure  6.  The  important  gap  size  is  then  the  apparent  trailing  edge 
height  h(t)  seen  by  an  inner  observer,  which  matches  with  an 
apparent  gap  width  as  x ”*  - •,  seen  by  an  observer  in  the  stern- 
edge  region,  i.  e. 


lim  fL(x,  t)  lim  f (x,  t) 

h(t)  ^x  — Xg-ot^  J=x~-oo  L ^ J . (A9) 

(x  = 0(1))  (x  - Xg  = 0(c)) 


Now  we  can  write  the  trailing -edge  solution  as 

x - x (t) 

0(x,y,t)  = 0g(t)  + h(t)U(t)  ®(  'l^t))  + (x-xg(t))Xg(t),  (A10) 

where  $(X,  Y)  is  the  canonical  steady-flow  potential  sketched  in  Figure 
6,  and  <t>^,  U are  to  be  determined.  Note  that  <b  possesses  a unit 
velocity  magnitude  as  X -*  + oo . We  assume  that  the  wake  for 
$(X,  Y)  is  bounded  by  a free  streamline  Y = F(X)  which  is  of  width 
o at  X = + co , relative  to  the  unit-width  channel  at  X = - oo  , and 
on  which  the  flow  generated  by  $ has  a unit  magnitude.  The  last 
term  in  (A10)  is  needed  to  bring  the  moving  body  surface  to  rest. 

The  problem  for  ®(X,  Y)  is  a classical  free -streamline  problem 
(c.  f.  Gilbarg,  1960),  and  hodograph  techniques  may  be  used  to  solve 
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for  $ , for  any  given  stern  geometry.  In  particular,  we  then  obtain 
the  coefficient  a , which  is  commonly  called  the  "coefficient  of 
contraction"  for  the  given  geometry.  Actual  contraction,  as  in  a 
nozzle,  applies  when  o < 1;  Figure  6 illustrates  an  "expansion" 
case  with  a > 1.  If  a rudder  is  present,  one  would  expect  o < 1 
when  the  rudder  is  to  port  and  o > I when  it  is  to  starboard.  The 
constant  C shown  in  Figure  6 can  also  be  computed,  but  plays  no 
essential  role,  to  leading  order  in  t . 

Verification  of  the  validity  of  (A10)  demands  only  that  <p  matches 

x - x. 

with  the  inner  flow  as  — pj *♦  - oc  , and  satisfies  the  unsteady 

wake-boundary  condition  (2.  10).  Clearly,  by  its  construction,  (t) 
is  identifiable  with  <p(x  (t),  t)  , and  since  * -*  a X as  X**  - oc  , 
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Substituting  for  U from  (All)  and  from  (A8),  we  have  (with 
<t>  = 0(xg.  t) , u s = <t>'  (xg,t))  , 


+ usxs  f — 2 («s 


1 • 2 
2 XS 


(A1 3) 


which  is  equivalent  to  (3.  14),  after  some  manipulation. 
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Figure  1. 
Figure  2. 

Figure  3. 

Figure  4. 


Figure  5. 

Figure  6. 


FIGURE  CAPTIONS 


Sketch  of  flow  and  geometry. 

Transient  wake  velocity  and  thickness  at  ^ = 3 , for 
an  airfoil  started  impulsively  at  t = o with  constant 
speed  U , plotted  against  x for  various  times  t . 
Transient  wake  as  in  Figure  2,  at  p - 3 . Note 
"shock"  (fw  - oo)  for  t/tQ  - 2,4. 

Scaled  height-versus -time  trajectories  for  fall  of  a 
plate  of  (scaled)  mass  m - 1,  for  various  (scaled) 
horizontal  speeds  U . Arrows  indicate  times  of 
impact  for  U = 0,1;  there  is  no  impact  if  U > 1.  64. 
Canonical  unit-velocity  flow  from  a unit -width  straight 
channel,  emerging  into  a half-plane,  as  a source  of 
strength  2 at  infinity. 

Canonical  free -streamline  flow  near  shaped  stern, 
or  stern  with  rudder.  Channel  flow  has  velocity  o 
and  width  1,  whereas  wake  at  infinity  has  velocity  1 
and  width  a , where  o is  the  contraction  coefficient 
of  the  effective  nozzle  so  formed. 
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